Astronomy & Astrophysics manuscript no. H4322 


February 2, 2008 


(DOI: will be inserted by hand later) 





The intrinsic colour dispersion in Type la supernovae 

S. Nobili^, A. Goobar^ R. Knop^, and P. Nugent^ 



cn 
o 
o 



> 
o 

(N 

O 
m 
o 

Oh 

6 



Department of Physics, Stockholm University, 

SCFAB, S-106 91 Stockholm, Sweden 

Department of Physics & Astronomy, 

Vanderbilt University, USA 

Lawrence Berkeley National Laboratory, 

1 Cyclotron Road, Berkeley, CA 94720, USA 



Abstract. The properties of low-redshift Type la supernovae are investigated using published multi-band optical broadband 
data from the Calan/Tololo and CfA surveys. The average time evolution of B-V, V - R, R-I,B-I and V - I, the intrinsic 
dispersion and time correlations are studied. This information is required to deduce the extinction of such explosions from the 
measured colours. We find that extinction corrections on individual SNe based on their colours up to 40 days past the fi-band 
lightcurve maximum are generaly limited to cr^|, >0.1, due to intrinsic variations, as far as it can be conservatively deduced 
with the current sample of data. However, we find that the V - R colour, especially at late times, is consistent with a negligible 
intrinsic spread, and may be the most accurate estimator for extinction. 
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1. Introduction 

In the last few years Type la supernovae (SNe la) have proved 
to be excellent distance estimators and have been successfully 
used t o investigate the f ate of the universe jPerlmutter et al.L 
Il999t iRiess eTall \l99?i h. Despite the broad use of these ob- 
jects by cosmologists, the cuiTent knowledge of the nature of 
SNe la is rather limited. Thus, there is a strong demand for fur- 
ther understanding to assess important issues about these ex- 
plosions. For cosmological implications, the main concerns are 
related to the possible evolution of the SN properties with red- 
shift. Moreover, critical tests for extinction along the line of 
sight based on supernova colours require good knowledge of 
the intrinsic properties of these objects. 

In this paper, a statistical study on 48 well observed nearby 
SNe la is caiTied out. In particular, the intrinsic dispersion in 
SN colours is investigated using published BVRI data. We also 
focus on the time correlation of intrinsic optical colours. This 
information is needed to address the possible host galaxy or 
intergalactic extincti on by dust of super n ovae used for cosm o- 
logical tests (see e.g. ' Riess et alJ ( l2000l) : iNobiU et alJ (l2003l) ). 
and also to probe for other exotic sources of dimin ing at high-z 
with differential extinction jMortsell et al. [Eo02h . 



2. The data set 

Published BVRI lightcurves of well observed nearby SNe la 
were analysed. The considered sample consists of 48 SNe la 
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from 2 di fferent sets, the Calan/Tololo data published by 
iHamuv e~al. ( 1996) , and the set in Riess et al, ( 1999 ), usually 
referred to as the CfA data. The list of SNe is given in Tabled 
along with the observed filter data available for each of them, 
th eir redshift and the jj-ba nd l ightcurve "stretch ", s, as defined 
in lPerlmutter et alJ ( Il997l) and lOoldhaber et alJ f2001 ). 
The selected samples include a broad variety of SNe, well dis- 
tributed in stre tch factor, s. This parameter, related to Ami5 
ilPhillir)stll99i . has been found to coiTelate with the supernova 
luminosity. Thus, a sample well distributed in stretch should 
imply a broad distribution in luminosity. Fig. ^ shows the dis- 
tribution of the Z?-band light curve stretch factor and SN red- 
shifts (0.003 < z < 0.12), for both samples. The timescale 
stretch para meter wa^determined from the lightcurve fits as in 
iGoldhaber et alJ ll200ll) using their B-band lightcurve template 
with a parabolic behavior for the earliest epoch after explosion. 

A'-corrections were applied to a ccount for the small cosmo- 
logical redshift as described in e.g. lGoobar et alJ( 2002). using 
the (s — 1) spectroscopic template of 'Nug ent et ah (120021) as 
a starting point. The results of this analysis were used to im- 
prove the spectral template of SN la's (Sect.|6|l and we iterated 
the analysis once re-calculating the /T-coiTections with the im- 
proved template. Note that even though the ^T-coiTections for 
the used data set are small, typically of the order of a few hun- 
dreds of a magnitude, for some of the more distant objects in 
the sample considered, they reach up to ~ 0.5 mag. 

The SNe Ught curves were corrected for both M ilky Way 

I 1 1 -1 

and host galaxy extinction as in Phillies et al. ( 1999) using the 

method first proposed bv .Lira L1995.) for estimating host galaxy 
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extinction using late epoch light curves. There is empirical ev- 
idence that the B - V colors of SNe la show extremely small 
scatter for the period between 30 and 90 days post B-band max- 
imum, despite any difference in the light curve shapes at earlier 
epochs. As the s et of SNe used in thi s article is a subset of the 
one analysed in IPhimns etal.l (11999'). the host galaxy extinc- 
tions listed in Table 2 of their paper were used. These were 
derived combining the late time B -V colour with information 

on the B - V and V - / at m aximum light. 

The lCardelli et al. I(ll989 ) relation, modified bv lO'DonneUl 

ill 994 . was used to compute the extinction in other colours 
given E{B - V). Spectral templates of Type la SNe were used 
to compute the evolution of the extinction with the supernova 
epoch. 

The extinction corrected lightcurves were further screened 
to exclude the most peculiar SNe, as the main emphasis of 
this work is to establish the properties of "normal" supernovae. 
Fig-Eshows the difference of the Z?-band lightcurve maximum, 
^max and the V-band light curve maximum V^ax plotted against 
the de cline rate parameter. Amis, as reported in Phillips et al] 
1I1999II) . A 3cr clipping rejection criteria was applied, iterating 
until no data points were further rejected. At least 5 SNe in the 
2 data sets de viate significantly f rom the expected linear rela- 
tion derived in lPhillips et al.l(ll99 9) from an independent set of 
20 non-reddened SNe (solid line in Fig.|2}: 

-0.07(+0.012) + 0.114(+0.037)(Ami5 - 1.1) 



Ta ble 1. List of SNe u sed f or the analysis 

rn. lRiess et al']lll998t) : (2\ \Uarmiv etaDJlw 



B, 



max ' max 



Vrr 



The outliers are SN 1993ae, SN 1995bd, SN 1996ai, 
and SN 1996bk for the CfA data, and SN 1992K for the 
Calan/Tololo set. Moreover SN 1995E and SN 1995ac have 
been excluded since they are similar to the peculiar SN 199 IT 
llBranchetal.[ll993l) .' 



Fig. 1. Distribution of redshift, z, (binsize=0.02), and stretch 
factor, s, (binsize= 0.08), for the analysed sample 



3. The average optical colours 

Figs. jSj- jS] show the time evolution of the extinction and K- 
con-ected B - V, V - R, R - I, B - I md V - I, colours. The 
plotted errors include the uncertainty on the host galaxy exinc- 
tion correction. The time axis (f') has been corrected by the SN 
redshift and the B-band stretch factor s, to account for the de- 
pendence of colours on the stretch. This rescaling of the time 



" excluded from the analysis presented here. 
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axis was found very effective in reducing the measured intrinsic 
dispersion (Sect.j^ 



' Including these 2 SNe in the analysis does not change the results 
significantely. 
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Fig. 2. Bmax-Vmax vs Amis for the 2 sets of data, k-corrected, 
corrected for Galactic exti nction and Host Gala xy extinction, 
compared with the result o flPhillips et alJ ( 1 19991) 



For each colour lightcurve we have spline-interpolated 
the weighted mean values computed in four days wide, non- 
overlapping bins. The aim of this procedure was to find a curve 
that describes the average time evolution of the colours. This 
parameterization will be referred as a "moder in the following 
discussion. Other methods were used to fit the data, as for ex- 
ample a least squares cubic spline fit. These result in curves that 
differ from the model typically by 0.01 magnitudes. The models 
show some systematic discrepancies between the two data sets, 
especially m R — I for which the C-T model is always redder 
than the CfA model. In Fig. |5]we investigate the differences 
of each of the models from the one built on both sets together. 
The differences are usually of the order of a few hundreds of 
a magnitude. The largest deviation was found for R - I in the 
Calan/Tololo set, resulting in a difference of about 0.2 mag at 
maximum with respect to the model built up using data from 
both sets. Note, however, that the statistics in the Calan/Tololo 
set for R - I colour, all along the evolution and in particular 
around the time of Bmax, is extremely poor, as shown in Fig.|5] 
Due to the smaller quoted observational error bars, the CfA set 
dominates the weighted average used to build the model out of 
both data sets. 

A source of uncertainty in this analysis is the ability of the 
observers to convert the instrumental magnitudes from the used 
filter-i-CCD system transmission into the standard BVRI sys- 
tem. The apparent systematic effects for R- l in Fig.|9lmay be 
indicative of this. 
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Fig. 2>. B -V for the two sets of data. The solid line represent 
alway the curve found for the current set, while the dashed line 
is the curve found for the same colour of the other set. 
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Fig. A.y - R for the two sets of data. The solid line represent 
alway the curve found for the current set, while the dashed line 
is the curve found for the same colour of the other set. 

In order to assess systematic effects in builing the colour 
models, we checked whether the host galaxy extinction was 
over-corrected. Thus, we comp ared the models with the data 
of those SNe that, according to Phillips et al. ('1999^, suffered 
no extinction from the host galaxy. The comparison, shown in 
Fig. exhibits no obvious deviations and we may conclude 
that the extinction corrected colours are consistent with the un- 
corrected sub-sample. 

A possible remaining dependence of colours on the stretch 
factor s was investigated. For data points in a 5 days broad 
bin around time of B -maximum and around day 15, a linear 
weighted fit was done on the residuals from the models not 
showing any further significant dependence on Amis or equiv- 
alently on the stretch parameter s. Note that this is not in con- 



4 



Nobili et al.: The intrinsic colour dispersion in Type la supernovae 



Calan/Tololo 




10 20 30 40 50 60 




10 20 30 40 50 60 



Fig. 5. R - I for the two sets of data. The soHd Une represent 
alway the curve found for the current set, while the dashed line 
is the curve found for the same colour of the other set. 



Calan/Tololo 




Fig. 6. B - I for the two sets of data. The solid line represent 
alway the curve found for the current set, while the dashed line 
is the curve found for the same colour of the other set. The filled 
circles represent SN 1993H, while the filled diamonds are data 
of SN 1992bc. Refer to Sect. for a discussion about these 2 
SNe. 

tradiction with Fig.|3 as we are considering (B- y),^^^^ instead 
of Bmax - VmiLx- Morcovcr, our analysis includes data points in a 
broad bin around day zero, and not a fitted lightcurve maximum 
as in Fig. 12 

4. Intrinsic colour dispersion 

The observed colour dispersion around the models is overall 
inconsistent with a statistical random distribution of the data 
with the reported measurement errors. The;^'^ value per degree 
of freedom computed for each model is typically between 3 
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Fig. 7. Comparison of the models with data of unreddened SNe 
only. 
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Fig. 8. V-I for the two sets of data. The solid line represent al- 
way the curve found for the current set, while the dashed line is 
the curve found for the same colour of the other set. The filled 
circles represent SN 1993H, while the filled diamonds are data 
of SN 1992bc. Refer to Sect.Qfor a discussion about these 2 
SNe. 



and 6, giving indication of a very poor fit of the data. Thus, 
we conclude that the analysed data supports the existence of an 
intrinsic colour dispersion of Type la supernovae ^. 

The intrinsic colour dispersion was computed on the resid- 
uals of each of the data sets from the corresponding mot/e/, built 



- Tiiis has been further investigated as described at the end of this 
section. 
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Fig. 9. Residuals of each of the models from the one built for 
all the data set. 

as explained above. LabeUng XYj any of the measured colours, 
B - V,V - R,R - I, B - I and V - I, the residual with respect 
to the model expectation is referred as Rxyi^ i e-j 



R 



XYi 



XYi - XY] 



model 



(2) 



Next, the time axis was binned into 5 days wide bins to com- 
pensate for the modest statistics. The bins are centered at days 
d - 5-n, with n ranging from to 8 . For each bin k, we compute 
the weighted average for the data XY, (see Appendix IXt 



(3) 



where Nk is the number of points in bin k and w, is the inverse of 
the uncertainty on the i-th measurement squared, l/aj, where 
o", include both measurement errors and uncertainty on the host 
galaxy extinction corrections. The weighted sample standard 
deviation was computed as the square root of the weighted sec- 
ond moment: 



.2 

XYi 



(4) 



The uncertainty on the expression in Eq. @ is given by the 

* 1 



square root of the variance, y[i^j,], computed as: 
VI^y] 



4m,, 



(5) 



where y[mw j is the variance of the weighted second moment. 
We can consider the result of Eq. @ as an estimate of the in- 
trinsic dispersion in each bin, for each XY colour: Axy - sty. 



This has been applied at each one of the sets separately and at 
the whole set of data. The results of the method for the whole 
set of data are given in Table|2| 
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Fig. 10. Comparison of the results for the intrinsic dispersion 
in each colour given in Table |2l for both data sets (triangles) 
and the results obtained keeping the 2 sets of data separated 
(squares and diamonds for the C-T and CfA respectively) 



Fig.^Jshows a comparison of the intrinsic dispersion com- 
puted for each of the sets separately and their combination. In 
most of the cases the differences exceed the statistical uncer- 
tainties, seemingly arising from systematic differences between 
the two data sets. The computed intrinsic dispersion for the CfA 
sample was found to be smaller than for the CT data set and the 
whole (combined) sample in most cases, again pointing at sys- 
tematic differences in the reported magnitudes. 

Note that the adopted method leads to an overestimation 
of the intrinsic colour dispersion due to the contribution from 
the measurement errors. However, the weighting procedure en- 
sures that the most accurate measurements dominate the re- 
sult. To asses the impact of the measurement accuracy, we 
run a Monte-Carlo simulation to generate a synthetic colour 
data set with a dispersion given by the measurement uncer- 
tainties alone, i.e. no intrinsic dispersion. Three hundred data 
sets, with the same distributions in epochs and formal error bars 
as the CfA and Calan/Tololo were simulated and the weighted 
standard deviation (and its error) were computed, according to 
Eqs. ^ and (|5jl. Note that the simulation, for simplicity, gen- 
erates gaussian distributed and completely uncorrelated data. 
The averages, 6, were used to disentangle the contribution of 
the intrinsic dispersion from the measurement errors. First, an 
hypothesis test was run to verify whether the simulated data 
and the measured data had the same dispersion; e.g. implying 
null intrinsic dispersion: 
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A level of significance a - 0.01 was set for rejecting the 
null hypothesis, ( Cowan, .1998). Only 10 cases were not re- 
jected, indicated by a in Table|2l For all the other cases, for 
which the Ho hypothesis was rejected, the intrinsic dispersion 
was computed as: 

A"'^ = VA2 - ^2 (6) 

and a lower limit on its value was set at a 99% confidence limit. 
The cases for which the null hypothesis was not rejected, were 
considered as compatible with no intrisic dispersion, and an 
upper limit on its value was set at a 99% confidence level. The 
corrected intrinsic dispersions are listed in the 5:th column of 
Table |2 together with upper and lower limits. We notice that 
the narrowest colour dispersion happens for V - R, especially 
at late times. At 25 days after Bmax and later, this colour is 
compatible with no intrinsic spread at all. Further, it should be 
noted that at day 35, all the colours but V- / are consistent with 
vanishing intrinsic dispersion. 

5. Correlation 

The correlation between optical colours at different epochs was 
also estimated. The property that was tested is whether a super- 
nova that is blue at a certain ep och for example, say at maxi- 
mum, stays blue at all epochs. In lRiess et al the authors 
argue that data measurements more than 3 days apart may be 
considered as uncorrelated estimators of colour. Our analysis 
does not support that assumption.-^ We find significant correla- 
tions for data points up to a month apart, as shown below. The 
method followed is essentially the one used to compute the in- 
trinsic dispersion. One can summarize the following steps: 

- Bin the data in time 

- Select only the SNe present in all the time bins 

- For each bin compute the weighted average of the measure- 
ments belonging to the same SN. 

- Compute the linear correlation coefficient between bins as 
in eauation lA-81 

- Test the correlation coefficient significance (ApDendix lA.2> 

We refer to appendix lA.2l for what follows. The correlation 
coefficient between different epochs h and k is: 

/j ^ 

rhk = , , _, (7) 

where the summation is on the i-th SN, which by construc- 
tion is present in all the bins. The uncertainty on ri,k has been 
computed converting it into the normally distributed variable z, 
as described in (Appendix IA.2> . The results, given in Table|3] 
show that the correlation is important and non-zero all along the 
time evolution. Fig. II II shows the weighted mean colours for 

^ For high-z objects this is even more questionable when some of 
the main sources of uncertainty is the subtraction of a common image 
of the host galaxy and the A^-corrections. 



the selected SNe in the bins centered at day 0, 15 and 30. Note 
that the supernovae selected are different in different colours, 
but, by construction, are the same in all the bins for each of the 
colour It appears that supernovae that deviate from the average 
colour at a certain epoch are likely to keep their colour excess 
all along the 30 days evolution considered here. 

6. Supernovae template 

iNugent et alJ ( l2002l) analysed the relation between colours of 
SNe la and /^-corrections. In particular they showed how K- 
corrections are mainly driven by the overall colour of the SN 
rather than by peculiarities of single features. Using spectra and 
colour light curves they give a recipe to build a SN la tem- 
plate, to be used for computing /T-corrections of ^ ~ 1 .0 SNe. 
As the sample used in our work is larger than the one consid- 
ered in their paper, we have used our results to improve the 
spectral template. Note that this affects only the BVRI mag- 
nitud es, whic h are the only bands treated here. Referring to 
Nuge nt et al ](|2002), we proceed as follows: 

- correct the BVRI of their list of UBVRIJHK magnitudes 
using our B - V, V - R and R - I models 

- correct the corresponding spectral templates with this 
colours'* 

- iterate the whole analysis described in this paper, using the 
newly created templates to compute the /T-corrections for 
both Calan/Tololo and Cf A data set 

- correct again the templates using the most recent version of 
the models 

Even though the template itself is modified considerably 
from the original one, iterating once results in small corrections 
(~ 0.001-0.005 mag) in the intrinsic dispersion in all colours at 
any epoch, and a few percent in the colour models. As the cor- 
rection obtained are small and well within the given uncertain- 
ties, we do not iterate further The results of this analysis should 
not he considered definitive and may change as SNe la and their 
colour evolutions are studied in even more detail. However, it 
can be considered a good estimate of the average SN la tem- 
plate (in BVRI), as it is constructed from a quite broad sample 
of SNe. The final corrected template is available upon request. 

7. Discussion 

Throughout our analysis all data points were treated as inde- 
pendent measurements. This is particularly important for the 
cases with significant host galaxy light underneath the super- 
nova where a single reference image was used, introducing a 
correlation between the data points not considered here. 

A first estimate of the intrinsic dispersion was calculated 
for each data set separately and the two data sets together A 
comparison of the results shows that the colours extracted from 
the CfA sample have smaller scatter indicating that the contri- 
bution from measurement errors is not negligible. Further, in 
the analysis the measurements are assumed to be Gaussian dis- 
tributed and the weighted standard deviation has been taken as 

The same code used bv lNugentetalJj2002h was used for this. 
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Table 3. Correlation coefficients between the different bins. The indicated xbin is the central value of each bin. The bin size is 7 
days for all the bins. The errors indicate the 1 cr confidence level for the computed coefficients. See text for details. 
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an estimate of the intrinsic dispersion. However, we noticed 
that 2 SNe, SN 1993H and SN 1992bc, seem to be rather de- 
viant in B - / and V -I for the C-T set, as shown in Figs. |6land 
|8l The effect of the "outliers" is particularly important aroud 
20 days after B-maximum. For comparison, we recalculated 
the intrinsic spread in B - / and V - I excluding these 2 SNe. 
The results are shown on Table |4] The agreement between the 
intrinsic dispersion between the data sets improves when these 
2 SNe are excluded. We emphasize that there are systematic 
differences between the Calan/Tololo and CfA data sets, and 
they might have been introduced while converting from the in- 
strumental system used to the standard BVRI system. 

An attempt of disentangling the intrinsic dispersion from 
the contribution of the uncertainties has been done resulting in 
upper limit values for the intrinsic dispersion in some cases. 
The most relevant is the case of V - R, that seems compati- 
ble with null intrinsic dispersion for most of the epochs at 99% 
C.L. Moreover this analysis brings out an important feature at 
day 35, when all the colours but V - I are consistent with zero 
intrinsic dispersion. This indicates that further studies will be 
needed to investigate this intriguing finding. When computing 
the correlation coefficients we considered data up to 30 days 
after B-band maximum, even though later epoch data are avail- 
able for several supernovae. However, the necessary condition 



of each SN being observed in all the time bins reduces the 
statistics if later epochs are introduced. This limitation is spe- 
cific to the sample used and can be overcome with the use of 
a more extensively observed sample, suc h as what will be pro- 
vide by SNfactory ('Alde ring et"all l2002V 

The intrinsic dispersion sets constrains on the ability to 
determine the host galaxy extinction,Av. This will depend on 
the colour, as the intrinsic dispersion is different for different 
colours. As an example we used the Cardelli et al. ( 1989) rela- 
tion at the effective wavelength for each bandpass to compute 
the expected uncertainty in the extinction, neglecting any de- 
pendence on the supernova phase. 





= 3.1 ■ 


A corr 
' ^BV 




= 6.2 


A COIT 
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Ay 


= 4.1 ■ 


A corr 
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= 1.4 ■ 
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"vi 



where Ry was assumed equal to 3.1. Table|5]shows the results 
of the equation |8l for the epochs for which the intrinsic disper- 
sion was calculated. The A^°" was used for this. The results 
indicate that, with the present knowledge, extinction by dust 
with Ry = 3.1 may only be determined to cr^^>0.10 with Type 
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Table 4. Results of the analysis done excluding SN 1992bc and SN 1993H (**) for B - I and V - I. xbin is the central value 
in days for each bin; XY^** is the weighted mean value and its uncertainty; A^/** is computed according to Eqs. and (|5|i; 
A™"** is the corrected intrinsic disperion and L.L is the lower limit at 99% C.L.; In the last column is the difference between the 
results given in Table |3and this analysis; the errors are computed as sum in quadrature. 
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la restframe optical data within the first 40 days after B-band 
lightcurve maximum, for the colours and epochs with non-zero 
intrinsic dispersion. However observations inV - R are prefer- 
able to other colours to set limit on the extinction of an ob- 
served supernova. To account for possible different extinction 
parameters Ry in the different s uper nova host galaxies (as n o- 
ticed in e.g. ('Ries s et aLl ll99(^ and (Krisciunas et al., 2000)), 
we repeated the analysis considering a Gaussian uncertainty 
= 1 on Rv, propagating this scatter on the host galaxy 
extinction corrections. This changes the values of the intrinsic 
dispersion, but typically within the quoted errors in Table|2] 

8. Conclusion 

A statistical analysis of colours of SNe la using a sample of 48 
nearby SNe was performed. Of this sample 7 SNe have been 
excluded based on their large extinction, or peculiar behav- 
ior of their colours at maximum. With the present knowledge 
and data quality we computed the average colour evolution for 
B - V,V - R,R - I, B - I and V - I with time and the derived 
intrinsic scatter. We find that the correlation of colours during 
the first 30 days after restframe B-band maximum is not negli- 
gible, i.e. arguing agai nst the assumptions made in the analysis 
of SN 1999Q (z^0.46: lRie'ss et alll200dl) where five measure- 
ments of restframe B - I along its lightcurve were treated as 
independent estimates of extinction. According to our findings, 
their limit on the presence of intergalactic grey dust must be re- 
vised. A reanalysis of the data from SN 1999Q is in preparation 
jNobili et al.. 2003) . With the data at hand, host galaxy extinc- 
tion corrections from restframe optical colours within the first 
30 days after maximum light are generaly limited to cr^,, >0. 1 
due to the intrinsic variation of Type la colours at those epochs. 



with the possible exception of extinction corrections derived 
from the rest-frame V - R colour The results of this analysis 
have been used to correct spectroscopic templates which are 
available upon request. 
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Appendix A: Appendix 

A.1. Algebraic moments 

In what follow s we refer to ICowanI ( Il998h and 
lKendall& Stuart I ( Il958h (see Volume 1, Chapter 10), 
our aim here is to give uniformity to the notations used. Given 
n independent observations of a variable, xi,X2, the r-th 
moment or algebraic moment is given by: 

OT, ^-V^ (A-1) 
n z—i 

(=1 

The expectation value, or mean value, and the variance of tUr 
are: 
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where jizr is the 2 ■ r-th moment. The IA-2l are exact formulae 
as long as one knows ^2r and Hr- However this is not always 
the case, and one has to use their estimators and from 
the sample itself. The variance of the standard deviation can be 
computed using standard error propagation: 



y[cr] = y[ v^] 



V[m2\ '«4 - '«2 



4mi 



Anrri') 



(A-3) 



In experimental situation each observation x, is often at- 
tached to a certain weight w,. Supposing that the weight them- 
selves are known without errors, one can define the following 
formula for the r-th weighted moments: 



(A-4) 



In order to evaluate the variance on the weighted standard de- 
viation we simply extended the IA-3l for the case of weighted 
moments, obtaining: 



4m„ 



An 



(A-5) 



where the m^.^ and m„., are the 4th and 2nd weighted moments 
respectively and n,.// = {Y!l=\ ^ifl T!l=\ 
A simple Monte Carlo simulation has been run to verify the 
accuracy of the approximated formula lA-51 

A.2. Correlation coefficient 

Given 2 random variables x and y the correlation coeflicient is 
defined as 



Pxy = 



(A-6) 



The unbiased estimator of the covariance Vxy is: 



— 1 " 

Vxy = - 'x){yi - y) 

n - 1 ^ 

1=1 



(A-7) 



SO that the estimator of the correlation coefficient will be: 



' xy 



Vxy 

SxSy 



ZU(xi-x)(yi-y) 



(A-8) 



For the IA-8l to be a good estimator of the correlation coef- 
ficient there are a few caveats to check. First condition is that 
the samples are randomly defined from the population, that is 
to make sure that the samples x and y are not selected in some 
ways that would operate as to increase or decrease the value of 
r. Even though one has random samples it is possible to com- 
pute the errors due to sampling. Commonly this is computed 
as (1 - r^)/ ■\/n. Unfortunately this is just an approximation. 
Moreover r's for successive samples are not distributed nor- 
mally unless n is large and the true value p is near zero. This 
yields to a distinction: 



if n > 30 

in order to know whether the value calculated for r is sig- 
nificantly different from zero, one can compute its standard 
error as: 



cr, = 1/ Vn- 1 



(A-9) 



If r/cr is greater than 2.58, one can conclude that the uni- 
verse value of r is likely to be greater than zero, 
if n < 30 
the variable: 



t - r- 



vr 



(A-IO) 



follows the t-distribution with dof=n-2. This can be used 
only for testing the hypothesis of zero correlation. 

R. A. Fisher developed a technique to overcome these diffi- 
culties. The variable r is transformed into another variable that 
is normally distributed. This is especially useful for high value 
of r, when none of the above test can be safely applied. The 
transformation to the variable z: 



z = iln(l+r)-iln(l-r) 



(A-11) 



allows some important semplifications. The distribution of z's 
for successive samples does not depend on the universe value p 
and the distribution of z for successive samples is so near to nor- 
mal that it can be treated as such without any loss of accuracy, 
see for example (Volume 1, Chapter 16 and Volume 2, Chapter 
26) of Kendall & Stuart ( 195^. Moreover the standard error 
for z is independent on its cr: 



cr, = 1/ Vn-3 

The way to proceed is then very simple: 



(A-12) 



- Compute r according to ( IA-8> 

- Transform r into z according to ( IA-1 1> 

- Compute cr, according to ( IA-12> 

- In order to state a 1 sigma confidence limit for r, transform 
the 2 values z + ctz back to r, using the inverse of HA-l U 

References 

Aguirre, A., 1999, ApJ, 512, L19 

Aguirre, A., 1999, ApJ, 525, 583 

Aguirre, A. & Haiman, Z., 2000, ApJ, 532, 28. 

Aldering, G. et al., 2002, Supernova Factory Webpage, 

http : //snf actory . Ibl . gov 
Branch, D., Fisher, A. & Nugent, P., 1993, AJ, 106, 2383 
Cardelli, J.A., Clayton, G.C. & Mathis, J.S., 1989, ApJ. 

345,245 

Cowan, G., 1998, Statistical data analysis, Oxford University 
Press 

Cramer, H., 1957, Mathematical methods of statistics, 

Princeton University Press, Seventh Printing 
Goldhaber, G., Groom, D.E., Kim, A. et al., 2001, ApJ, 558, 

359 



10 



Nobili et al.: The intrinsic colour dispersion in Type la supernovae 



Goobai-, A., Mortsell, E., Amanullah, R. et al., 2002 A&A, 392, 
757 

Hamuy, M., Phillips, M.M., Suntzeff, N.B. et al., 1996,AJ, 1 12, 
2408 

Kendall, M.G. & Stuart, A., 1958, The advanced theory of 
statistics, Charles Griffin & Company Limited, London. 

Krisciunas, K., Hastings, N.C., Loomis, K. et al., 2000, ApJ, 
539, 658 

Lira, P., 1995, Master Thesis, Univ. Chile 

Mortsell, E.,Bergstrom, L & Goobar, A., 2002, Phys. Rev. D, 

66, 047702 
Nobili, S. et al., 2003, in prep. 

Nugent, P, Kim, A., & Perlmutter, S., 2002, PASP 1 14, 803 
O'Donnell, J.E., 1994, ApJ, 422, 1580 

Perlmutter, S., Gabi, S., Goldhaber, G. et al., 1997, ApJ, 
483,565P 

Perlmutter S., Aldering, G., Goldhaber, G. et al., 1999, ApJ, 

517,565P 
Phillips, M.M., 1993, ApJ,413L.105P 

Phillips, M.M., Lira, P., SuntzefF, N.B., Schommer, R.A., 
Hamuy, M.& Maza, J., 1999, ApJ. 118, 1766-1776 

Riess, A.G., Press, W.H., & Kirshner, R.P, 1996, ApJ, 473, 
588R 

Riess, A.G., Filippenko, A.V., ChalHs, P et al., 1998, AJ, 116, 
1009. 

Riess, A.G., Kirshner, R.P,Schmidt, B.P et al., 1999, ApJ. 
117,707-724 

Riess, A.G., Filippenko, A.V., Liu, M.C.et al., 2000, ApJ, 536, 
62 



Table 2. Results of the analysis of all SNe. First column: cen- 
tral value in days for each time bin; A^^ is the number of points 
for each bin; XF* is the weighted mean value and its 1-sigma 
uncertainty; Axy is the intrinsic dispersion computed according 
to Eqs. (0} and (|5j; A™" is the corrected intrinsic dispersion, 
computed as in Eq. (|6}, and in the last column is the estimated 
lower limit at 99% C.L. (see text); 

" compatible with null intrinsic dispersion at 99% C.L.; an up- 
per limit is given instead of the corrected intrinsic dispersion. 
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Table 5. cr^,, represent the constrain on the extinction Ay that it 
is possible to compute for a given intrinsic dispersion. The A™^ 
and it L.L given in Table|3has been used for these estimates. 
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Fig. 11. From top to bottom and left to right. Residuals from the models for the SNe selected for each colour, in bin centered at 
day 0,15 and 30. The band defined by the solid lines correspond to the intrinsic dispersion found in at the same epoch. 



